Adiabatic state preparation of interacting two-level systems 
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We consider performing adiabatic rapid passage (ARP) using frequency-swept driving pulses to 
excite a collection of interacting two- level systems. Such a model arises in a wide range of many body 
quantum systems, such as cavity QED or quantum dots, where a nonlinear component couples to 
light. We analyze the one-dimensional case using the Jordan- Wigner transformation, as well as the 
mean field limit where the system is described by a Lipkin-Meshkov-Glick Hamiltonian. These limits 
provide complementary insights into the behavior of interacting systems under ARP, suggesting our 
results are generally applicable. We demonstrate that ARP can be used for state preparation in the 
presence of interactions, and identify the dependence of the required pulse shapes on the interaction 
strength. In general interactions increase the pulse bandwidth required for successful state transfer, 
introducing new restrictions on the pulse forms required. 



Precise control of quantum mechanical systems is a 
sought after feature for appUcations in quantum informa- 
tion and investigations of quantum dynamics. Discrete 
atomic-Uke systems, or qubits, can be excited by using 
external field pulses that induce Rabi oscillations, with 
the final state determined by the intensity and duration 
of the pulse. However, this method is sensitive to fluc- 
tuations in the driving field, transition energy and other 
sources of disorder Ij. An alternative approach, which 
is robust against such variations, is the use of frequency- 
swept ( "chirped" ) pulses to perform adiabatic rapid pas- 
sage (ARP). In this method, the frequency of the driving 
field is swept through the transition to be excited, imple- 
menting the Landau-Zener process for adiabatic passage 
[21 13] . Provided the gap induced by the applied field is 
large compared with the sweep rate the process is adia- 
batic, and the wavefunction is transferred from the initial 
ground state to the target state with high probability. 

ARP is a well-established technique in nuclear mag- 
netic resonance, where chirped radio frequency pulses 
are used to manipulate nuclear spins [4]. More recently, 
there have been a number of investigations into using 
ARP with optical pulses to control excitons in quantum 
dots [TJ HHT] , including the creation of entangled states 
[SlllOj. This has coincided with growing interest in pro- 
ducing many-body systems with strong light-matter in- 
teractions, such as coupled photon cavities or polaritonic 
systems TP. A protocol such as ARP that allows ro- 
bust control of the quantum state in these systems would 
enable the investigation of quantum dynamics in highly 
non-equilibrium regimes |12H15j . 

In established examples of ARP the interactions are 
weak on the scale of the level splittings generated by 
the ARP pulse, and hence the former can be straight- 
forwardly neglected. The aim of this paper is to demon- 
strate how ARP may be extended to strongly-interacting 
regimes where this is not the case. We consider a model 
of interacting two-level systems which, by comparison to 



the case of a single two- level system [T1, allows the ef- 
fect of interactions to be identified. We show that ARP 
remains an effective approach in the interacting case, 
provided the pulse bandwidth is sufficient to span the 
spectrum of the collective modes generated by the in- 
teractions. Although our model is relatively simple, our 
results are relevant across a wide range of systems, in- 
cluding cavity QED systems mj , quantum dots [101 dU j 
superconducting qubits [17l [18] and doped impurities in 
semiconductors |19j . 

The model we consider consists of a set of N interact- 
ing two-level systems driven by an external field (in the 
rotating wave approximation): 
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where <Tj are the Pauli matrices for the two-level system 
«, and af = (af ± iaf)/2. In this form the two states 
af — ±1 are understood to correspond to the presence 
or absence of an excitation of the «th two-level system, 
e.g., of an exciton in a particular state of a particular 
quantum dot. We will also refer to the collective pseu- 
dospin S ~ Tlii'^i/'^^ whose z-component is related to 
the total excitation or occupation n = S'^ -f N/2. fi{t) 
is the coherent external pulse used to perform ARP, and 
Jij is the interaction between systems i and j. At this 
stage we assume that the energy of the excitation E > Q 
is the same for all transitions, and neglect interactions of 
the form cF^a^. This model could be realized in precisely 
engineered cavity [TT] or circuit QED [T71 [TB] systems. 
Furthermore, a less idealized model of this form can be 
used to describe many realizations of interacting qubits, 
such as coupled quantum dots [TOl [16]. These systems 
often exhibit disorder in the energies E and interaction 
strengths Jij , but the robustness of ARP means the gen- 
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eral understanding we obtain of the effect of interactions 
is applicable. 

Decomposing the driving field into amplitude and fre- 
quency fi(t) = Qilt) exp(j j uj{t')dt'), and eliminating the 
instantaneous frequency from the driving term using a 
unitary transformation, the Hamiltonian becomes: 
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For the discussion in this paper, we consider a Gaussian, 
linearly chirped pulse with uniform amplitude, 
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where g parametrizes the pulse amplitude, r is the tem- 
poral width of the pulse and a is the linear chirp. We 
discuss the pulse and system parameters in terms of the 
dimensionless combinations gr, Jt and olt^ = !)• For 
a = 0, Eq. (|3]) becomes a Rabi pulse centered at fre- 
quency iJ, with a pulse area proportional to gT. 

In the non-interacting case where J = 0, the use of 
ARP to transfer a two-level system from the ground state 
cr^ = — 1 to the excited state is well understood [1 . For 
(7 = 0, the energies of the two levels cross when E + 
a;(t) = 0. The presence of the field g 7^ produces 
an anticrossing and the adiabatic state smoothly varies 
from the initial ground state to the excited state. When 
the pulse amplitude is time independent, g{t^ = 5, the 
model reduces to the canonical Landau-Zener problem 
[2l|3j. The probability of remaining in the adiabatic state 
(and so being transferred from the initial ground state to 
the excited state) is 1 — exp(— 27rg^/a), so that the final 
population is always increased by reducing the chirp a, 
increasing the adiabaticity of the process. In the case 
of ARP, using pulses of finite duration, g(t) is no longer 
constant. Thus, in order for adiabatic passage to occur, 
the two levels of the system must be coupled together 
long enough that the character of the eigenstates changes 
sufficiently slowly. This introduces the requirement that 
a ^ p]. In the limit a — 0, the system undergoes 
Rabi oscillations rather than ARP. 

In order to understand how this process generalizes to 
the interacting case, we first examine a one-dimensional 
chain with nearest neighbor interaction Jij = JSi^i+i. In 
this case, the energy levels for g = can be determined 
using the Jordan- Wigner transformation erf — 2clci — 1, 
a~ = exp(i7r ^^^j cjcj)ci = T^q where Ci are fermionic 
operators j^Oj- After also performing a Fourier transfor- 
mation the Hamiltonian, Eq. ([2]), becomes: 
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where A'' is the number of sites and k = —Tr+2TTTn/N with 
m integer. The Jordan- Wigner transformation has pre- 
viously been used to describe Landau-Zener transitions 
for anisotropic spin chains in a changing magnetic field 
pn [22]. The Landau-Zener transitions in that model 
result from the anisotropy, which affects the subspaces 
spanned by fermion operators of a given |A;| indepen- 
dently. In contrast, the spatial dependence of the non- 
linear Ti term in our model leads to terms in the Hamil- 
tonian that couple fermion states with different |fc|. 

In the Jordan- Wigner representation the different en- 
ergy eigenstates correspond to different occupations of 
the fermion states. The completely empty (spin down) 
state corresponds to the vacuum with no fermions |0). 
Likewise, the completely occupied (spin up) state corre- 
sponds to the case with all fermion states filled, Ylk "^l 
For large \t\ the first term in Eq. Q dominates and the 
eigenstates, shown in Fig. [T] for A^ — 4, are split into 
N + 1 bands labelled by the total number of fermions, 
which physically corresponds to the excited-state popu- 
lation of the two-level systems n. 

In the non-interacting case, the energy levels are in- 
dependent of k and all the states in the nth band have 
energy —nat. The presence of interactions, J 7^ 0, lifts 
the complete degeneracy of states within each band as 
shown in Fig. [l] The separate states correspond to the 
different allowed values of the total spin for a given S^. 
In order to prepare a fully occupied state, the quantum 
state must then be transferred via multiple level crossings 
from the n = to 71 = Af bands 

The splitting of the level crossings that allows adia- 
batic state transfer is introduced by the external pulse 
field g. In the general case there will be some variation 
in the driving field between the two-level systems, which 
makes the form of the coupling term in the fermionic 
representation complicated to determine due to the non- 
locality of the Jordan- Wigner string T^. However, for the 
uniform driving we consider [gi{t) = g{t)], the coupling 
to the field in the non-transformed Hamiltonian, Eq. ([2| , 
can be rewritten git)Y,^{a+ + a~) = g{t){S+ + S^), 
and the transitions therefore conserve S'^ . If the system 
starts in the ground state, it will thus always be in an 
eigenstate of with the maximal eigenvalue S{S + 1) 
where S — N/2. There is one state in each band with 
this value of S^, and the transitions between these states 
have matrix elements gA/(]v'^^ny[(7r+T72y^l72l. The 
corresponding states in each band are the most symmet- 
rical states, which for J > (J < 0) have the lowest 
(highest) energies, see Fig. [T] (Fig. [T] inset) pS] . 

The field term in the Hamiltonian only changes the 
number of fermions, n, by ±1. An anticrossing between 
non- adjacent bands can be induced by higher order vir- 
tual transitions. For example, a n — ^ (n + 2)th band 
transition is possible via an intermediate (n -I- l)th band 
state. These higher-order interactions are suppressed in 
the mean-field limit discussed below [231 US] i but do play 




FIG. 1. (Color online) Time-dependent eigenenergies of the 
one-dimensional chain, Eq. Q, with g — 0, J = ^ and = 4 
sites. Different colors indicate bands corresponding to dif- 
ferent values for the occupation n of the two-level systems. 
Bold lines show the states in each band with maximum , 
which are coupled together when the pump is spatially uni- 
form. Dashed lines are other eigenstates of the system, which 
have different values of . The vertical dotted lines show 
the separation in time between the crossing of the n — band 
with the n — 1 (A) and n = N (B) states. Inset: As main 
figure but with J = — |. Note that the order of crossings of 
the coupled (bold) levels has reversed. 



a role in ARP if the connectivity is small. 

The Jordan- Wigner transformation is only usefully ap- 
pHcable for the special case of nearest-neighbor hopping 
in one dimension. In higher dimensions, an alternative 
approach is to use the mean field approximation, which 
is exact in the limit N ^ oo, Jij = J/N"^. We show 
in Fig. [2] results for the final occupation obtained using a 
spatially uniform pulse, Eq. ([s]), calculated by solving the 
Heisenberg equations of motion using the mean field re- 
placement J2ijJij'^t^J = J2iJcsicrt i'^t^) + h.c). In 
this approximation the Hamiltonian, Eq. ([2|, can be 
rewritten in terms of the collective spin operators as the 
Lipkin-Meshkov-Glick Hamiltonian: 
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where we have used = {S' f + {S+S-+S-S+)/2and 
dropped terms which do not affect the dynamics. This 
Hamiltonian, with a time independent g{t) = g, has been 
used to describe Landau-Zener tunneling for a bosonic 
Josephson junction [24 l I26ti30] . 

Figure [2] consists of a fan of non-zero occupation 
bounded by large regions of essentially zero occupation. 
Within the fan, as in the normal Landau-Zener prob- 
lem, increasing a decreases the final occupation as the 




FIG. 2. Mean-field calculation of the average excitation of a 
set of interacting two- level systems, Eq. ([2|, driven from its 
ground state by the chirped pulse, Eq. (|3| with gr = 3. Jr is 
the dimensionless interaction strength, and ar the dimension- 
less chirp. The regions where no excitations are created are a 
result of the finite duration of the pulse r: as the chirp a is re- 
duced, the level crossings of Fig. [l] no longer occur within the 
pulse and so adiabatic transfer is not possible. The boundary 
of this region is approximately a oc J. 



increased velocity of the level crossing reduces the adia- 
baticity of the transition. 

For a fixed value of a, the occupation within the fan 
increases (decreases) for J > (J < 0). This varia- 
tion corresponds to the changing relative positions of the 
level crossings, visible in Fig. [T] for the one-dimensional 
chain. In the absence of interactions, J = 0, the states 
within each band are degenerate and so all levels cross 
simultaneously a.t t — 0. As the interaction strength J is 
increased, the level crossings separate in time. As each 
level degeneracy becomes more isolated, the size of the 
anticrossing caused by g increases. When J > and 
a > the crossings occur in "ascending order", i.e. the 
n ^ n + 1 crossing occurs before the n + 1 — n + 2 
crossing. The increase in splitting due to the isolation 
of crossings then improves the efficiency of transfer to 
the occupied state, producing the increase of the occu- 
pation shown in Fig. [2j If J < 0, however, then for 
positive chirp the crossings occur in the "wrong order" 
(see Fig. [T] inset), so that it becomes more difficult for 
the system state to transfer via a series of transitions 
through adjacent bands, suppressing the probability of 
full occupation. In the mean field limit, this ordering 
leads to the formation of a swallowtail in the energy level 
evolution which causes a breakdown of adiabaticity and 
a corresponding reduction in the occupation [231 . 

The large regions of zero occupation that define the fan 
are a result of the time dependence of the optical pulse 
g{t) used to perform ARP. In order for the state to be 
transferred at a level crossing, the anticrossings caused 
by the field g{t) must be large enough compared to the 
level velocity a to make the process adiabatic. Because 
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ARP uses pulses of a finite duration there is only a lim- 
ited window during which g{t) meets this criterion. With 
no interaction, J = 0, all crossings occur simultaneously 
and so may all occur within the window. In the presence 
of interactions, eventually the time of the first level cross- 
ing J/a ^ T, and it will be pushed out of the pulse. Ne- 
glecting the change in anticrossing size discussed above 
and the effect of higher order, virtual transitions, the 
crossover to non-adiabatic behavior occurs along a line 
J (X a, which is approximately what is seen in Fig. [2] 
Below this line, the behavior is no longer adiabatic and 
the system undergoes more complicated dynamics, re- 
ducing to non-linear Rabi oscillations for a = 0. Thus a 
pulse used to perform ARP in a large system must have 
a sufficient duration t that it includes the entire region 
of level crossings separated by the interaction. 

In conclusion, we have shown the consequences of inter- 
system interaction on using ARP to fully occupy a set 
of two level systems. The interaction lifts the degener- 
acy of the eigenstates of the Hamiltonian, Eq. Q, and 
causes the level crossings at which state transfer occurs 
to separate in time. As in the bosonic Josephson junc- 
tion PH[?7H5^ . the isolation of each degeneracy increases 
the effective splitting so that adiabatic transfer can be 
achieved for larger chirps than in the non-interacting sys- 
tem. However, the separation of level crossings also in- 
troduces the additional condition for ARP that the pulse 
duration (or chirp) should be large enough to include all 
the necessary crossings, increasing with J. Physically, 
this occurs because the interactions broaden the spec- 
trum into a set of collective modes forming a path from 
the ground to final states, and the pulse must cover this 
spectrum for the state preparation to be effective. 

Although in this paper we have focused on the ideal 
case of uniform E, J and coupling gi, our results apply 
more generally. Fluctuations in E and J will change the 
energies and character of the intermediate states so that 
they are not delocalized across the system. However, for 
\t\ — 7> oo, the highest and lowest states remain the empty 
or full states, so our results will still apply. Variation in 
gi changes the size of splittings at a level crossing but, 
with the exception of fine-tuned cases, anticrossings will 
still form, allowing adiabatic transfer. 

As the model discussed in this paper represents lim- 
its of more complicated many-body systems including 
the Bose-Hubbard, Dicke or Jaynes-Cummings-Hubbard 
models [HI \17\ |3T] , these results can be used as a basis 
for understanding the behavior of ARP in these models. 
It could then be used as a robust method of preparing 
far-from-equilibrium states in those systems for use in 
quantum information contexts or as equivalents of the 
quantum quenches performed in ultracold atomic gases. 
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